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ABSTRACT 

Aims. We aim to reveal environmental dependence of infrared luminosity functions (IR LFs) of galaxies at z~0.8 using the AKARI 
satellite. AKARFs wide field of view and unique mid-IR filters help us to construct restframe 8/im LFs directly without relying on 
SED models. 

Methods. We construct restframe 8^m IR LFs in the cluster region RXJ1716.4+6708 at z=0.81, and compare them with a blank 
field using the AKARI North Ecliptic Pole deep field data at the same redshift. AKARFs wide field of view (10' x 10') is suitable to 
investigate wide range of galaxy environments. AKARFs 15/im filter is advantageous here since it directly probes restframe 8/im at 
z ~0.8, without relying on a large extrapolation based on a SED fit, which was the largest uncertainty in previous work. 
Results. We have found that cluster IR LFs at restframe 8/im have a factor of 2.4 smaller L* and a steeper faint-end slope than that 
of the field. Confirming this trend, we also found that faint-end slopes of the cluster LFs becomes flatter and flatter with decreasing 
local galaxy density. These changes in LFs cannot be explained by a simple infall of field galaxy population into a cluster. Physics 
that can preferentially suppress IR luminous galaxies in high density regions is required to explain the observed results. 

Key words, galaxies: evolution, galaxies:interactions, galaxies:starburst, galaxies:peculiar, galaxies:formation 



1. Introduction 

It has been observed that galaxy properties change as a function 
of galaxy environment; the morphology-density relation reports 
that fraction o f elliptical galaxies is larger at higher galaxy den- 
sity dGoto et a l., 2003) ; the star fomiation rate (SFR) is higher in 
lower galaxy density dGomez et all 120031: iTanaka et al.L |2004) 
. However, despite accumulating observational evidence, we 
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still do not fully understand the underlying physics governing 
environmental-dependent evolution of galaxies. 

Infrared (IR) emission of galaxies is an important 
probe of galaxy activity since at higher redshift, a sig- 
nificant fract ion of star fo rm ation is obscured by dust 
(Takeuchi, B uat. & Burgarellal 120051; Goto et al., 201_d. 
Although there exist low-z c luster studies ( Bai et all 120061 ; 
IShim et all l2010t iTran et all l2010h . not much attention has 
been paid to the infrared properties of high redshift cluster 
galaxies, mainly due to the lack of sensitivity in previous IR 
satellites such as ISO and IRAS. Superb sensitivity of recently 
launched Spitzer and AKARI satellites can revolutionize the 
infrared view of environmental dependence of galaxy evolution. 
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Fig.l. Restframe 8^im LFs of cluster RXJ1716.4+6708 at 
z=0.81 in the squares, and those of the AKARI NEP deep 
field in the triangles. For RXJ1716. 4+6708, only photometric 
and spectroscopic cluster member galaxies are used. For the 
NEP deep field, galaxies with photo-z/specz in the range of 
0.65 < z < 0.9 are used. The dot-dashed lines are 8/im LFs 
of RXJ1716. 4+6708, but scaled down for easier comparison. 
The thin dotted lines are the best-fit double power laws. Vertical 
arrows show the 5cr flux limits of deep/shallow regions of the 
cluster (red) and the NEP deep field (blue) in terms of Lsf_mi at 
z=0.81. 



In this work, we compare restframe 8/im LFs between clus- 
ter and field regions at z=0.8 using data from the AKARL 
Monochromatic restframe 8/im luminosity (is/im) is important 
since it is kno wn to correlate well with the total IR luminosity 
jBabbedge et al.. 2006; Huang et al., 2007), and hence, with the 
SFR of galaxies (Kennicutt. 119981) . This is especially true for 
star-forming galaxies because the rest-frame 8/im flux are dom- 
inated by prominent PAH features s uch as at 6.2, 7.7 and 8.6 fj,m 
(iDesert. B oulanger. & Puget'.'l99(/). 

Important advantages brought by the AKARI are as follows: 
(i) At z=0.8, AKARI's 15/im filter (L15) covers the redshifted 
restframe 8/im, thus we can estimate 8/fm LFs without using 
a large extrapolation based on SED models, which were the 
largest uncertainty in previous work, (ii) Large field of view of 
the AKARI's mid-IR camera (IRC, 10' x 10') allows us to study 
wider area including cluster outskirts, where impo rtant evolu- 
tionary mechanisms are su ggested to be at work too to et all 
I2004t iKodama et all llool . For example, passive spira l galax- 
ies h ave been observed in such an environment (Got o et al.L 
I2003h . Unless otherwise stat ed, we adopt a cosm ology with 
{h,nm,^A) = (0.7,0.3,0.7) dKomatsu et a"n.l2009l) . 

2. Data & Analysis 

2. 1. LFs of cluster RXJ1 71 6.4+6708 

The AKARI is a Japanese infrared satellite (iMurakami et al.L 
I2007h . which has continuous filter coverage in the mid 
IR wavelengths (iV2, 7V3, iV4, S7, S9W, 511, L15, L18W and 
L24). The AKARI has observed a massive galaxy cluster. 



Fig. 2. Restframe S/zm LFs of cluster RXJ17 16.4+6708 at 
z=0.81, divided according to the local galaxy density {T,^th)- 
The stars, circles and squares are for galaxies with logYi^th > 2, 
1.6 < logT,r-,th < 2, and logT,r-,th < 1-6, respectively. 

RXJ1716.4+6708, in iV3, 5*7 and L15 (iKovama et al.L l2008h . 
RXJ1716.4+6708 is at z=0.81 and has a = 1522tji5kin s"!, 
Lx,„, = 13.86 ± 1.04 X lO'^^ erg s-\ kT = 6.8^1° keV. Mass 
estimate from weak lensing and X-ray are 3.7+ 1.3 x 10^"^ Mg 
and 4.35+0.83 x IO^'^Mq, respectively (see iKovama etatl 
|2007 » for references). 

An important advantage of the AKARI observation is L15 
filter, which corresponds to the restframe 8/im at z=0.81. With 
15 (3) pointings, L15 reaches 66.5 (96.5)/iJy in deep (shal- 
low) regions at 5a. Here flux is measured in 11" aperture, 
and coverted to total flux using AKARI's IRC correction table 
(2009.5. ifl Cluster studies with the Spitzer are often performed 
in 24/xm and thus needed a large extrapolation to estimate ei- 
ther Lsf_mi or total infrared luminosity (Ltir, 8 — lOOO/iw). 
Note that we do not claim the Lg^m is a better indicator of 
the total IR luminosity than other indica tors jBrandl et al.Ll2006t 
ICalzetti et al 1 l2007t iRieke et al.L |2009|) . but it is important that 
the AKARI can meausure redshifted 8/iTO flux directly in one of 
the filters. 

Thanks to the AKARI's wide field of view (10' x 10'), the 
total area coverage around the cluster is 200 arcmin^, which 
cover larger area than previous cluster studies with the Spitzer, 
allowing us to study IR sources in the outskirts, whe re important 
galaxy evolution takes place (e.g.. lGoto et al.Ll2003l) . Previously, 
^oyama eVa.1. (2008) reported a high fraction of L15 sources in 
the intermediate density region in the cluster, suggesting a pres- 
ence of environmental effect in the intermediate density environ- 
ment. 

This same region was imaged with Suprime- Cam in VRi'z' 
and h as a good photometric redshift estimate jKovama et al.L 
I2007h . Used in this work are 54 L15-detected galaxies which are 
well identified with optical sources with 0.76 < Zphoto < 0.83. 

With the L15 filter covering the restframe 8/im, we simply 
convert the observed flux to 8/im monochromatic luminosity 

' http://www.ir.isas.jaxa.jp/ASTRO- 
F/Observation/DataReduction/IRC/ApertureCorrection_09050 1 .html 
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Table 1. Best double power-law fit parameters for LFs 
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Sample 




0*(Mpc~^dex-') 


a 


/3 


NEP Deep field 
RXJ17 16.4+6708 


6.1±0.5 X 10'" 
2.5±0.1 X 10'° 


0.0010±0.0003 
0.74±0.04 


1.1±0.3 
2.6±0.1 


5.7±1.2 
5.5±0.4 



(^S/im) using a standard cosmology. Completeness was mea- 
sured by distributing artificial point sources with varying flux 
within the field and by examining what fraction of them was re- 
covered as a function of input flux. Since we have deeper cover- 
age at the center of the cluster, the completeness was measured 
separately in the central deep region and the outer regions of 
the fie ld. More detail of the method is described in IWada et al.l 
(l2008h . 

Once the flux is converted to luminosity and completeness is 
taken into account, it is straight forward to construct Lspm LFs, 
which we show in the squares in FigU] Errors of the LFs are 
assumed to follow Poisson distribution. Here, we take an angular 
distance of the most distant source from the cluster center as 
a cluster radius {Rmax = 6.2Mpc). We assumed ^t^R^^^ as 
the volume of the cluster to obtain galaxy density (</)). This is 
only one of many ways to define a cluster volume, and thus, a 
caution must be ta ken to compare abs olute values of our LFs to 
other work such as Shim et al. ('2010'). This cluster is elongated 
in angular direction ( Koyama et al., 2007) , and thus, the volume 
might not be spherical. Yet, comparison of the shape of the LFs 
is valid. 

2.2. LFs in the AKARI NEP Deep field 

Our field LFs are based on the AKARI NEP Deep field 
data. The AKARI performed deep imaging in the North 
Ecliptic Pole Region (NEP) fr om 2-24/xm, with 4 po i ntings 
in each field ov er 0.4 deg^ (Mat suhara et all l2006l 120071: 
IWada etaLl |2008|) . The 5 a sensitivity in the AKARI IR filters 
(N2, m, Ni, S7, S9W, Sll, Lib, LISW and L24) are 14.2 , 
11.0, 8.0, 48, 58, 71, 117, 121 and 275/iJy (IWada et alLlIOOl . 
Flux is measured in 3 pix radius aperture (=7"), then corrected 
to total flux. 

A subregion of the NE P-Deep field (0.25 deg^) has ancillar y 
data from Subaru BV Ri' z' \lma.\ si al.Ll2007l;IWada et aUl2008l) . 
CFHT u'(Serjeant et al. in prep.), KPN02m/FLAMINGOs J 
and Ksdmm et al., 2007), GALEX FUV and NUV (Malkan 
et al. in prep.). For the optical identi fication of MIR sources, 
we ad opt the likelihood ratio method ([Sutherland & Saundersl 
119921) . Using these data, we estimate photometric re dshiftof L15 
detected sources in the region with the LePhare dllbert et al.L 
[2006; Arnouts et al., 2007). The measured errors on the photo-z 
against 293 spec-z galaxies from Keck/DEIMOS (Takagi et al. 
in prep.) are ^^^=0.036 at z < 0.8. We have excluded those 
sources better fit with QSO templates from the LFs. 

To construct field LFs, we have selected Lib sources at 
0.65 < Zphotoz < 0.9. There remained 289 IR galaxies with 
a median redshift of 0.76. Lib flux is converted to Lsf^m us- 
ing the photometric redshift of each galaxy. LFs are com- 

Siuted using the l/Vmax method. We used the SED templates 
Lagache. Dole. & Pugetl 120031) for /c-corrections to obtain the 
maximum observable redshift from the fl ux limit. Completenes s 
of the Lib detection is corrected using iPearson et all ( l2009bl) . 
This correction is 25% at maximum, since we only use the sam- 
ple where the completeness is greater than 80%. 

The resulting field LFs are shown in the dotted line and tri- 
angles in FiglT] Errors of the LFs are computed using a 1000 



Monte Caro simulation with varying z and flux within their er- 
rors. These estimated errors are added to the Poisson errors in 
each LF bin in quadrature. 

We performed a detaild com parison of restframe 8/im 
LFs to those in the literature in Goto et al. (2010). Briefly, 
there is an oder of dif ference between Caputi et al. ( 2007) and 
iBabbedge et al.l (l2006l) . reflecting difficulty in estimating Ls/™ 
dominated by PAH em issions using Spitzer 24/im flux. O ur field 
8/im LF lies between ICaputi et~al] (|2b07) and Babbedg e et al.l 
(2006). Compared with these work, we have directly observed 
restframe 8/im using the AKARI Lib filter, eliminating the un- 
certaintly in flux conversion based on SED m odels. More detail s 
and evolution of field IR LFs are described in lGoto et al.l ( l2010l) . 

3. Results & Discussion 

3. 1. 8fim IR LFs 

In FiglT] we show restframe 8/im LFs of cluster 
RXJ1716.4+6708 in the squares, and LFs of the field re- 
gion in the triangles. First of all, cluster LFs have by a factor 
of ^^700 higher density than the field LFs, reflecting the fact 
the galaxy clusters is indeed high density regions in terms of 
infrared sources. 

Next, to compare the shape of the LFs, we normalized the 
cluster LF to match the field LFs at the faintest end, and show 
in the dash-dotted line. In contrast to the field LFs, which show 
flattening of the slope at logLg/im < 10. SL©, the cluster LF 
maintains the steep slope in the range of 10. Oi© < logLg^m < 
10.6^0. The difference is significant, considering the size of er- 
rors on each LF. 

We fit a double-power law to both cluster and field LFs using 
the following formulae. 



$(L)dL/L* = $*|^— j dL/L*, {L<L*) (1) 
/ r \ 1-/3 

= $*f — j dL/L*, {L>L*) (2) 

Free parameters are: L* (characteristic luminosity, Lq), cf)* 
(normalization, Mpc^'^), a and /3 (faint and bright end slopes), 
respectively. The best fit values for field and cluster LFs are sum- 
marised in Table[T]and shown in the dotted lines in FiglT] 

The bright-end slopes are not very different, but L* of the 
cluster LF is smaller than the field by a factor of 2.4, and the 
faint-end tail of cluster LF is steeper than that of field LF. 

To further examine the difference at the faint end of the 
LFs, we divide the cluster LF us i ng the local galaxy density 
(Est/i) measured bv lKovama et al.l (|2008|) . This density is based 
on the distance to the 5th nearest neighbor in the transverse di- 
rection using all the optical photo-z members, and thus, is a 
surface galaxy density. We separate LFs using similar criteria, 
log'S^th > 2 (dense), 1.6 < logY^^th < 2 (intermediate), and 
logY.Qth < 1-6 (sparse), then plot LFs of each region in the 
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stars, circles, and squares in Fig|2] A fraction of the total vol- 
ume of the cluster is assigned to each density group in inversely 

3 /2 

proportional to the sum of Eg^,^ of each group. 

Interestingly, the faint-end slope becomes flatter and flatter 
with decreasing local galaxy density. This result is consistent 
with our comparison with the field in FiglT] In fact, the lowest 
density LF (squares) has a flat faint-end tail similar to that of the 
field LF. Since these LFs are based on the same data, changes in 
the faint-end slope are not likely due to the errors in complete- 
ness correction nor calibration problems. The completeness of 
the deep and shallow regions of the cluster are measured sep- 
arately. The changes in the slope is much larger than the maxi- 
mum completeness correction of 25%. We also checked the clus- 
ter LFs as a function of cluster centric radius, to find no signifi- 
cant difference, perhaps due to the elongated morphology of this 
cluster At the same time, assuming the same cluster volume, 
Fig|2] shows that a possible contamination from the field galax- 
ies to cluster LFs is only ^^0.1% in the dense region and 
even in the sparse region. 

It is interesting that not just the change in the scale of the 
LFs, but there is a change in the L* and the faint-end slope (a) 
of the LFs, resulting in the deficit in the 1O.2L0 < logig^m < 
10. 8^0 for cluster LFs. One might imagine a change just in L* 
might explain the difference in FiglT] However, in Fig|2] there 
clearly is a change in the slope as a function of H^th- 

However, interpretation is rather complicated; a shape of LF 
would not change if field galaxies infall into cluster uniformly 
without changing their star-formation activity. Although in clus- 
ter enviro nment, a fraction of M IR luminous galaxies is smaller 
than field (iKovama et al.L l2008l) . uniform and instant quenching 
of star-formation activity of field galaxies can only shift a LF, 
but cannot account for a change in L* and a of the LFs. 

Two important findings in this work are; (i) L* is smaller 
in the cluster (ii) the faint-end slopes become steeper toward 
higher-density regions. To explain these changes in LFs, IR- 
luminous galaxies need to be preferentially reduced, with a rela- 
tive increase of IR-faint galaxies. However, an environmental- 
driven physical process such as the ram-pressure stripping or 
galaxy-merging would quench star-formation not only in mas- 
sive galaxies but in less massive galaxies as well, and thus is not 
able to explain the observed changes in LFs. 

On the other hand, it has been frequently observed that more 
massive galaxies formed earlier in the Universe. This downsiz- 
ing scenario also depends on the environment, in the sense that 
galaxies with same mass are more evolved in higher density en- 
vironments than galaxis in less d ense environments (IGoto et al.L 
120051; iTanaka et al.1 12005L |2008[) . Statistically, a good correla- 
tion h as been found between Ltir and stellar mass (Elbaz et al., 
12001) . Our finding of the relative lack of IR-luminous galaxies 
in the cluster environment may be consistent with the downsiz- 
ing scenario, where higher density regions have more evolved 
galaxies and lacks massive star-forming galaxies. In contrast, 
in lower density regions more massive galaxies are still star- 
forming. However, since the data we have shown is in IR lumi- 
nosity, to conclude on this, we need good stellar mass estimate 
based on deeper near-IR data. 

Although a specific mechanism is unclear, the steep faint- 
end could also result from the enhanced star-formation in less 
massive galaxies. In the above scenario, massive galaxies have 
already ceased their star-formation in the cluster, but less mas- 
sive galaxies are still forming stars. These less massive galaxies 
may stop star-formation soon to join the faint-end of the red- 
sequence (iKovama et al.Ll2007l) . 
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Fig.3. Total infrared LFs of cluster RXJ1716.4+6708 atz=0.81 
in the solid line, and those of the AKARI NEP deep fi eld in the 
dashed line. Overplotted are the LFs of MS 1054 from lBai et al.l 
(l2007h. 



3.2. Total IR LFs 

To compare the Ls/xm LF in FiglUto those in the literature, we 
ne ed to convert Lsum to Ltir- We use the the following relation 
bv lCaputi etaTI (120071) ; 



L 



TIR 



= 1.91 X (i.L,,^,Ji-"6(±55%) 



(3) 



This is better tuned for a similar luminosity ra nge used here 
than the original relation by Bavouzet et al.l (l2008l) . The conver- 
sion, however, has be en the largest sourc e of errors in estimat- 
ing Ltir from Ls^.m. lCaputi et al.l (l2007l) report 55% of disper- 
sion around the relation. It should be kept in mind that the rest- 
frame 8/im is sensitive to the star-formation activity, but at the 
same time, it is where the SED models have strongest discrep- 
anci es due to the compl icated PAH em ission lines (see Fig. 12 of 
Caputi.etaLl[2007; Got o et aUlMoh . 

The estimated Ltir can be, then, converted to SFR using the 
following relat ion for a Salpeter IMF, (j) (m) cx m^^^^ between 
0.1 - IOOMq (lKennicuttlll998h . 



SFR{MQyr-^) = 1.72 x lir^^ Ltir{Lc, 



(4) 



In Fig|3] we show the Ltir LFs. Symbols are the same as 
Figin In the top axis, we show corresponding SFR. Overplotted 
asterisks are cluste r LF o f MS1054 at z=0.83 with x2 larger 
mass by iBai et al.l (l2007l) . which show good agree ment with 
our LF s of RXJ17 16.4+6708 in the squares. Note that lBai et all 
(l2007h covered only the central region of MS 1054 due to the 
smaller field of view of the Spitzer The shape of their LF looks 
more similar to our LFs in the highest density bin in Fig|2] A 
shift in scale is perhaps due to difference in esimating cluster 
volumes. 

A major difference of our work to that o f iBaiet al.l(l2007l) IS 
that they were not able to compare in detail on the shape of the 
LFs between field and cluster regions, due mainly to a smaller 
field coverage and larger errors on LFs. They had to fix the faint- 
end slope with a local value. The largest source of errors is in 
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converting Spitzer 24/im flux into 8/im. Both cluster and field 
LFs of this work use L15 filter, which measures restframe 8/im 
flux directly, eliminating the largest source of errors. In addition, 
both cluster and filed LFs are measured with an essentially same 
methodology, allowing us a fair comparison of LFs. 

4. Summary 

We constructed restframe 8/im LFs of a massive galaxy cluster 
(RXJ1716.4+6708) and a rarefied field region (the NEP deep 
field) at z^O.8 using essentially the same method and data from 
the AKARl telescope. AKARJ's 15/im filter nicely covers rest- 
frame 8/im at z~0.8, and thus we do not need a large interpola- 
tion based on SED models. AKARI's wide field of view allows 
us to investigate variety of cluster environments with 2 orders of 
difference in local galaxy density. 

We found that L* of the cluster 8/im LF is smaller than the 
field by a factor of 2.4, and the faint-end tail of cluster IR LFs 
become steeper and steeper with increasing local galaxy density. 
This difference cannot be explained by a simple infall of field 
galaxies into a cluster. Physics that preferentially supresses IR 
luminous galaxes in higer density regions is needed to explain 
the observed results. 
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